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Abstract 

A potential flow model is derived for a large flap-type oscillating wave energy 
converter in the open ocean. Application of the Green's integral theorem in 
the fluid domain yields a hypersingular integral equation for the jump in 
potential across the flap. Solution is found via a series expansion in terms 
of the Chebyshev polynomials of the second kind and even order. Several 
relationships are then derived between the hydrodynamic parameters of the 
system. Comparison is made between the behaviour of the converter in the 
open ocean and in a channel. The degree of accuracy of wave tank exper- 
iments aiming at reproducing the performance of the device in the open 
ocean is quantified. Parametric analysis of the system is then undertaken. 
It is shown that increasing the flap width has the beneficial effect of broad- 
ening the bandwidth of the capture factor curve. This phenomenon can be 
exploited in random seas to achieve high levels of efficiency. 

1 Introduction 

Recent analysis [1] has shown that the oscillating wave surge converter (OWSC) 
is an effective contender in the challenge to extract energy from water waves. 
In its simplest form, the OWSC is a bottom-hinged buoyant flap used to ex- 
tract wave power in the nearshore environment, where water particles move 
essentially in surge. 



Traditional research in wave energy has mostly focused on two-dimensional 
(2D) devices [2J [3J see for example], on three-dimensional floating bodies in 
deep water [2] and on buoys of characteristic width much smaller than the in- 
cident wavelength, known as "point absorbers" (PA) jU [5] . Regrettably, the 
potential of the OWSC has been somewhat disregarded in the past. Given 
the incident wave wavelength A and the flap width w, traditional theories 
are indeed applicable to the OWSC in the 2D approximation tc>A and in 
the PA limit w <C A [3 El see for example]. However, there is still a gap 
of knowledge in the region w = 0(A), which is the most likely for practical 
applications [lj. Lately, Falnes and Hals |7J proposed to treat wave energy 
converters (WECs) for which w = 0(A) as "quasi-point absorbers" (QPA). 
However, such definition is somehow ambiguous: a precise mathematical for- 
mulation is needed to determine the hydrodynamic behaviour of such WECs 
in detail. The mathematical models of Renzi and Dias [H [9] partially achieve 
this goal. They describe the resonant behaviour of a single OWSC in a chan- 
nel and of a periodic array of OWSCs. Within the framework of a potential 
flow theory, Renzi and Dias [HI E] show that the OWSCs can attain high 
levels of capture factor, with peaks occurring at resonance of the transverse 
short-crested waves of the channel/array system. Despite the latest advance- 
ments, however, the problem of determining the hydrodynamic behaviour of 
a single OWSC in the open ocean is still unsolved. 

In §2] a semi-analytical model for a single OWSC in the open ocean is 
derived, assuming the fluid inviscid and incompressible, the flow irrotational 
and the perturbation time-harmonic. The equation of motion of the flap is 
solved in the frequency domain and the hydrodynamic characteristics of the 
system are determined. The potential field at large distance from the flap 
is then analysed, to obtain the asymptotic expression of the radiated and 
diffracted waves in the far field. Several expressions are derived, which relate 
the hydrodynamic characteristics to the wave amplitude of the radiated and 
diffracted waves at infinity. In £J3] the model is validated against available 
numerical results. Then discussion is made on the main parameters of the 
system, for a layout inspired by the design of a commercial OWSC. In §H 
comparison is made between the behaviour of the OWSC in the open ocean 
and in a channel. This allows to quantify the degree of approximation of 
experimental models in wave tank, aiming at reproducing the performance 
of the converter in the open ocean. Finally, in £J5] parametric analysis of the 
OWSC in the open ocean is undertaken to investigate the influence of the 
flap width on the performance of the device. 
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Figure 1: Geometry of the system in physical variables; (a) plan view, (b) 
section. 

2 Mathematical model 
2.1 Governing equations 

Referring to figure [Jja, b), consider an OWSC in an open ocean of constant 
depth h'\ primes denote physical quantities. The OWSC is modelled as a 
rectangular flap of width w' and thickness 2a' . Under the action of monochro- 
matic incident waves of amplitude A'j and period T", the flap oscillates about 
a hinge placed upon a bottom foundation of height d . Since the device is 
designed to work in the nearshore, where the incident wave fronts are al- 
most parallel to the coastline, normal incidence is also assumed. Let 6'(t') 
be the pitching amplitude of the device, positive if anticlockwise and let t' 
denote time. Consider a Cartesian system of reference 0'(x', y', z'), with the 
X ~£IX1S opposite to the direction of propagation of the incoming waves, the 
y' axis along the width of the device and the z' axis pointing upwards from 
the still water level. The origin O' is in the middle of the device at the still 
water level (see again figure [TJx, b). Let us assume that the fluid is inviscid 
and incompressible and the flow is irrotational. Hence there exists a velocity 
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potential $>'(x', y', z', t') which satisfies the Laplace equation 

V 2 &(x',y',z',t') = Q (1) 

in the fluid domain; V'/' = (f' x /, f' y >, f x i) and subscripts with commas denote 
differentiation with respect to the relevant variable. Let us also assume 
that the incident wave amplitude is small as compared to the width of the 
device, which is the reference length scale of the problem: A'j/w' < 1. As a 
consequence, the angular rotation of the flap induced by the incoming waves 
is small and the behaviour of the system is linear [8]. On the free surface, 
the linearised kinematic-dynamic boundary condition reads 

^+9^ = 0, *' = 0, (2) 

with g the acceleration due to gravity. Also, absence of flux at the bottom 
of the ocean requires 

$; 2 , = o, z' = -ti, (3) 

while absence of normal flow through the flap yields 

& x , = -e' t ,(t'){z' + ti -c')H(z' + ti -c'), x' = ±0, \y'\ < w'/2, (4) 

where the thin-body approximation has been applied [H [101 see] . In expres- 
sion P|, the Heaviside step function H is used to model the absence of flow 
through the bottom foundation. Finally, the wave field generated by the in- 
teraction betwe en the inc oming waves and the flap must be outgoing at large 
distance r' = \/ x' 2 + y' 2 from the body. Let us now introduce the following 
non-dimensional variables 

(x,y,z,r) = (x',y',z',r')/w', t = Jg~J^>t ', $ = (^i'A')- 1 ®', = (w'/A')9' 

(5) 

and constants 

(h,c) = (h',c')/w', A I = A' I /A', (6) 

where A' — 1 m is the amplitude scale of the incident wave. By making 
use of §5§ and ([H]), the governing equation ([T]) and the associated boundary 
conditions ©-((l]) become, respectively, 

V 2 $ = (7) 

in the fluid domain, and 

®,tt + $ , z = 0, z = 0, (8) 
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and finally 



$ z = 0, z= -h, 



(9) 



& JX = -0,t(t)(z + h-c)H(z + h-c), x = ±0,\y\<l/2. (10) 
2.2 Solution 

Let us assume that the flap performs sinusoidal oscillations of frequency 
lo = u'^/w'/g = 2tt/T, for which the time dependence of the system variables 
can be separated as 

0(f) =a{6e- <wt } ) $(z,y,*,*) = K{(0 / + D + i? )e-^}. (11) 
In the latter, 

^ = iAl COshk ( Z + h ) c -ikx ^ 

to cosh kh 

is the spatial potential of the plane incident wave [6], while <p D (x,y,z) and 
(f> R (x,y,z) are, respectively, the potential of the diffracted and the radiated 
waves. In (fl2"l) . is the wavenumber, solution of the well-known dispersion 
relationship u> 2 = k tanh kh. By substituting the decomposition ([IT]) for the 
velocity potential into the governing equations (!7|)- (|T0|) . the latter become, 
respectively, 

V 2 (R,D) = (13) 

in the fluid domain, the shorthand notation 0( H,Z) ) denoting either potential, 
and 

0(f^=O, z = -h, (15) 

and finally 




x = ±0, \y\ < 1/2. (16) 



In the latter, V = iooQ is the angular velocity of the flap, depending on the 
complex amplitude of rotation (see fTTj) . still unknown. The expressions in 
f|T6|) are motivated by the fact that <p R is the solution of the problem in the 
absence of incoming waves, while (p D is the solution of the problem where the 
flap is held fixed in incoming waves [01 see]. Finally, (j)^ R,D ^ must be outgoing 
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as r — > oo. Following the procedure of [8], let us factor out the vertical 
dependence of the governing system (|T5|) - (fTB|) by setting 



ra=0 

where 



s,y,*) = ^P D) (*,y)S„(*), (17) 



a/2 cosh K n (2 + 

ZnW = 7I 2 • u2 ru72 » n = 0,l,... (18) 

( /i + a; smh K n fiJ 

are the vertical eigenmodes of the flat-bottom water wave problem [5J see], 
orthogonal along the water column: 

/ Z n (z)Z m (z)dz = S nm , n,meN, (19) 
J-h 

being S nm the Kronecker symbol. In ([TBI . k = fc, while k„ = zA; n are the 
solutions of the dispersion relation 

u 2 = —k n taiak n h, n — 1,2,... (20) 

By substituting the decomposition (TTTj) into (|T3|) and (|T6|) and using the 
orthogonality relation (fl9|) . we obtain the Helmholtz equation 



(V 2 + <p™>) = (21) 
and the boundary conditions 



x = ±0, |y| < 1/2, (22) 



for the unknown two-dimensional velocity potentials (p^' D \ where n = 0, 1, . . . 
is omitted from now on, for the sake of brevity. In expression (I22l) . 



Jn - ~ : : . o 77172 



y2 [k„(/i — c) sinh k„/j + cosh /t n c — cosh K n h\ 
(h + u~ 2 sinh 2 K n h) 

and 



k (h + 00 2 sinh 2 k„/i) 1//2 , 

dn = " ] 5 n (24) 

V2o; cosh K n h 
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are real constants. Finally, the y?n R ' D ' ) must be outgoing disturbances as r — > 
oo. The boundary- value problem fl2Tl) - fl22l) set for (pn' D ^ is similar in form to 
that obtained by Renzi and Dias [8] for an OWSC in a straight channel. In 
the latter, however, waves are reflected by the lateral walls and can propagate 
only along the channel axis. Following the general method of [8], the solution 
to (12~T1) - (I2~2"]) is found in [A] by applying an integral-equation technique based 
on the Green's theorem. Using an appropriate decomposition of the 3D 
Green function into a singular part and an analytic remainder allows to 
solve the plane boundary- value problem for tpn and to express the spatial 
potentials ^ R ' D) (PU) in a semi-analytical form (see \M for details). The 
radiation potential is given by 

— iV °° P 
(p R (x,y,z) = —r-^2K n xZ n (z)^2a {2p ) n 



X 



n=0 p=0 

1 



fay/a? + (y -u/2)*) 
^/x 2 + (y-u/2f 



(1 - u^U 2p {u) \ „ - dum 

j-i 

where HrP is the outgoing Hankel function of first kind and order n and 
the U 2 p denote the Chebyshev polynomial of the second kind and even order 
2p, p = 0, 1, . . . , P G N. Finally, the «(2 P )o are the complex solutions of the 
system of equations (152]) . ensuring that the first equation of ( !T6|) is satisfied. 
The «(2p)n are obtained with a numerical collocation scheme, therefore the 
solution ( 12 5 p is semi-analytical. Analogously, the spatial diffraction potential 
is given by 

p 

<f) D (x,y,z) = —^-kxZ (z)^2/3 {2 p)o 

(k^ + (y-u/2)A 
(1 - u 2 f 2 U 2p (u) \ } du. (26) 

In the latter, the /3( 2 p)n are the complex solutions of the system of equations 
(|62|) . ensuring that the second expression of ([IE]) is satisfied. Let us now 
analyse the equation of motion of the flap. 
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2.3 Flap motion 

In expression fl25|) the angular velocity of the flap, V = iuQ, is still unknown 
and must be determined by studying the motion of the body. The equation 
of motion of the flap in the frequency domain is that of a damped harmonic 
oscillator [5J, i.e. 

(-cu 2 I + C- iuupto) 6 = 7, (27) 

where oo is the frequency of oscillation, / = I'/(pw' 5 ) is the moment of 
inertia of the flap, C = C'/(pgw' A ) is the flap buoyancy torque and v vto = 
u ' P to/ (p w>5 \fgjw') is the power take-off (PTO) coefficient; p is water density. 
In expression (127)1 . 

ro yi/2 

7 = ibj \ \ A<p(y,z){z + h-c)dydz (28) 

J -h+c J -1/2 

is the complex hydrodynamic torque on the flap, where 

(-0, y, z) + fl (-O, y, z)} - U D (+0, y, z) + cf> R (+0, y, z)] (29) 



is the jump in potential between the two sides of the flap (see lAl) . By using 
( |29l) . ( 1571) and ( 1611 in ( |28|) . substituting the latter in ( 1271) and performing 
some algebra, the equation of motion of the flap reads 

[-u 2 (I + p) + C-iu(v + v pto )]e = F, (30) 

where 



is the exciting torque, 



-^Ajf (3 00 (31) 



U=o j 

is the added torque due to inertia and finally 

^=^-/oS{«oo} (33) 

is the radiation damping. In the latter, only the propagating mode n = 
is present, because higher order eigenmodes correspond to evanescent waves 
that do not radiate from the flap (see §2.41) . In the following, approximated 
expressions of the radiation and diffraction potentials will be obtained at 
large distance from the flap. 
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2.4 Behaviour in the far field 



First, let us consider the radiation potential f l25|) and rewrite it in polar 
coordinates: (p R (x,y,z) = <p R (r, 7,2), where 

(x,y) = r(cos7,shi7) (34) 

(see figure [TJi). For an observer far away from the OWSC, the flap behaves 
dipole-mode radiator, so that 

<j) R (r,"/,z) ~(f) R (r,0,z) cos j (35) 

as r — > 00. Now perform the change of variable (I34p into (I2l)j) and consider 
the inner integral of R (r, 7 = 0, z), i.e. 



H? K n rJl + £ 



I np = I (1 - u 2 ) 1 ' 2 ^^) V , J -du. (36) 

u 2 

4r 2 



For n > 0, then n n = ik n and the Hankel function inside the integral ( )36l) 
transforms into the modified Bessel function of the second kind K\ (TTJ see 
§8.407], which decays exponentially for large r. Therefore the terms n > do 
not give any significant contribution to the radiated wave field (|35|) at large 
distance from the flap. Physically, such terms correspond to the evanescent 
modes of the system, which are trapped near the source [SI see for example] . 
Now consider n — 0. At large r the Hankel function in Iq p ( 1361) can be 
approximated as [TQl see] 



^ 1) |/-\/l + 4^ 



v? 



u 2 



exp [ ikr\/l + — - «3vr/4 ) . (37) 



Note that the argument of the square roots in the latter expression is 1 + 
0(r~ 2 ). As r 00, the terms 0(r~ 2 ) in ( 137|) and in (|36|) can be neglected, 
thus yielding the sought approximated expression 



1 



ttI / 2 ifer 



2 r V 7rfcr 



^ /4 5 P o, (38) 
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Figure 2: Approximation of the integral I 00 (!36|) (solid blue line) with ex- 
pression ( 13 8 p (dashed black line) for large r. Parameters are: (a) k = 2, (b) 
k = 4. The geometry of the system is detailed in table [U 

as r 4 oo, where the orthogonality relation of [11, §7.343] has also been 
used. Numerical investigation for a typical geometry of the system reveals 
that fl38|) converges to the exact value fl36|) . as shown in figure [2]^ a, b). Now, 
expressing (l25l) in polar coordinates ( |34l) . applying the dipole relation ( |35l) 
and using the approximated expression ( |38l) for the inner integral, yields the 
expression of the radiation potential in the far field 



as r — oo. In the latter expression, 

.p. . — in , co cosh kh /A . 

^ (7 ) = _, cos7aM __^ (40 ) 

is the angular variation of the radially spreading wave [6l see], for unit ro- 
tational velocity of the flap. Analogous reasoning applied to the diffraction 
potential <p D ( 1261) yields 



a; coshfc/i V 7r«r 

at large distance from the plate, where 

.n, \ —in , n w cosh kh , 

7 = ^fccos 7 /3 o T79 42 

8v/2 (/i + ^sinh 2 ^) 172 
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is the angular variation of the diffracted wave, for unit amplitude of the 
incident wave. In the following section we shall derive some useful relations 
for the hydrodynamic coefficients of the system. 

2.5 Derivation of relations for the OWSC in the open 
ocean 

In this section we obtain some relations between the hydrodynamic coeffi- 
cients of the system, based on the asymptotic analysis of §2.41 Some of these 
relations are well-known in the literature for 3D bodies of general shape, 
while some others are new and provide insightful information on the hydro- 
dynamics of the OWSC. 

2.5.1 Relation between F and .4.^ (3D Haskind relation) 

Consider the exciting torque F fl3"T|) and rewrite it as 

F = ^pUjaoodo (43) 

by virtue of ( |65|) . where do is still given by ( 1241 . Then isolating «oo from ( HOI ) 
with 7 = 0, substituting it into (j4"3"l) and performing some simple algebra 
yields 

F = ^A I C g A R (0), (44) 

which is the well-known Haskind relation for a 3D floating body [6j see for ex- 
ample]. The latter relates the exciting torque to the wave amplitude radiated 
by the body in the direction opposite to the incident wave. 

2.5.2 Relation between F and A D 

Consider again the exciting torque F fl3Tl) . Now isolate the term /3oo from 
(j42p . substitute it into fl3Tl) and perform the straightforward algebra to obtain 

^ 4 , ,d, Aa.nh.kh / cosh kc — cosh kh\ , . 

= k ^ k ( ft ~ c + ksinhkh )■ («) 

The latter shows that the exciting torque is also proportional to the amplitude 
of the wave scattered in the opposite direction to that of the incident wave. 
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2.5.3 Relation between v and A 



Consider now the radiation damping coefficient v (133]) . By isolating aoo in 
f )40|) with 7 = and substituting it into fl33|) . the latter gives 

4<vx r ,p M i tanh&Zi /, cosh kc — cosh . . 

V = ft{A"lO)}—^-(h-c + kamhkh j, (46) 

which relates the radiation damping to the amplitude of the radiated wave 
in the forward x direction. 



2.5.4 Relation between u and F 

Finally let us inspect expression (1461) . By taking the real part of the Haskind 
relation (144}) . isolate 3? {.4.^(0)} and substitute it into (l46j) . so that the latter 
becomes 

tanh kh / cosh kc — cosh „ f F 1 , , 

v = -scr{ h ~ c+ kdnhkh ) X US- (47) 

which can be regarded as a particular form of the general Haskind relation 
( jH| tailored to the OWSC. Expressions (j4"4"]) - (l4"Tl) have been used to check 
the numerical evaluations of [A] Given the general equation l.h.s = r.h.s, the 
relative error 

|l.h.s. — r.h.s. | 
|r.h.s.| 

is defined to assess the accuracy of computations. For a typical OWSC 
configuration (see tabled]), taking maxn = 2 and maxp = 5 in (IB2"]) allows 
to obtain a maximum relative error e = O(10~ 15 ) for expressions (|44p - ()47|) . 
In the following section we shall analyse the performance of the device. 



2.6 Wave energy extraction 

Consider the equation of motion of the flap (]3"0l . The latter describes a 
damped harmonic oscillator in the frequency domain. As a consequence, the 
average generated power over a period is 

\F\ 2 

P = -TT^ r> (48) 



12 



provided the PTO system is designed such that 



Vo - v V-^\ s, (49) 

which corresponds to the optimum PTO damping [51 see for example]. To 
assess the performance of the OWSC, the capture factor 

P \F\ 2 
° F = \A]C g = 2A]C g {u pto + u) (50) 

is defined as the ratio between the power output of the device and the power 
of the incident wave per unit width of the device, where 

u ( 2kh \ 

is the group velocity of the incident wave in non-dimensional variables. Fur- 
thermore, if the flap is designed such that its buoyancy torque C and inertia 
torque / satisfy the relation 

u 2 (I + fi) = C, (52) 

i.e. the flap is tuned to resonance with the incident waves, v vto = v from ( H9l) 
and the capture factor (|50|) attains its maximum value 



~ AAjCgis m{A D (0)Y K } 

where (|45p and (|47[) have also been used. Note that the maximum capture 
factor (153|) depends on the amplitude of the diffracted wave in the direction 
opposite to that of the incident wave and can be determined directly by solv- 
ing the diffraction problem for (p D (|26|) . In the following, the mathematical 
model presented in this section will be validated against available numerical 
results. Then discussion will be made on the hydrodynamic coefficients of a 
typical OWSC configuration. Afterwards, the behaviour of a single OWSC 
in the open ocean will be compared to that of an OWSC in a channel, which 
reproduces a common layout for experimental studies in wave tank. Finally, 
parametric analysis will be performed to investigate the influence of the flap 
width on the performance of the OWSC in the open ocean. 



13 



Flap width (m) 


Foundation height (m) 


w 1 


d 


18 


1.5 


Ocean depth (m) 


Amplitude incident wave (m) 


h! 




10.9 


0.3 



Table 1: Geometry of the system and parameters of the incident wave field 
chosen for comparison with numerical data. 

3 Discussion 
3.1 Validation 

In this section we shall compare the results of the mathematical model of $2] 
with available numerical data. The geometry of the system chosen for com- 
parison is reported in table [U Such dimensions are inspired by the design 
of the Oyster 1™ wave energy converter developed by Aquamarine Power 
Ltd [12, APL, see]. Values of C and I have been obtained by private com- 
munication with APL. Note that, since A'j/w' ~ 0.02 <C 1, the regime of 
the system is linear. In figure |3] (a, b) comparison is made between the semi- 
analytical results of §2] and the numerical data of van't Hoff [13] for a single 
OWSC in the open ocean. The latter have been obtained with the software 
package WAMIT [13] for the linear analysis of fluid-structure interaction. 
In [13] the flap was modelled as a rectangular, 18 m wide box standing on 
a triangular prism pointing downwards, hinged on a 1.5 m tall platform in 
a 10.9 m deep ocean. Figures G^a) and E^fe) show, respectively, the values 
of the added inertia torque fi' = ppw' 5 and the magnitude of the exciting 
torque \F'\ = \F\pgA'w' 3 against the period of the incident waves for the 
semi-analytical model of §|2] and the numerical model of [13] . The agreement 
between both models is very good for the exciting torque and satisfactory 
for the added inertia torque. For the latter, minor discrepancies (see figure 
I3h) are likely due to the difference in shape between the flaps used in the 
models. Let us now analyse the performance of the device. 
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Figure 3: Comparison of the semi-analytical model of $2] (solid lines) with 
the numerical model of van't Hoff [13] (dots), for several thicknesses of the 
plate 2a'. (a) added inertia torque (132]) ; (6) exciting torque ( 131]) . Values are 
in physical variables. 



3.2 Performance analysis 

Figure HJ^a, b) shows the magnitude of the rotation |0'| = (Aj/u>')|0| and the 
capture factor Cp for the configuration of table [Q While the angular dis- 
placement of the flap increases monotonically with the incident wave period 
T', the capture factor instead is larger in short periods. The largest capture 
factor is max {Cp} ~ 0.73 at T 1 ~ 5.7s and is of the same order of magnitude 
as that predicted by Whittaker and Folley pQ for a similar OWSC configu- 
ration. For real devices, more complex flap shapes and foundation systems 
could help to achieve even larger figures [15J. Note that the behaviour of 
the capture factor of figure H^o) resembles very closely that of the exciting 
torque shown in figure [3)^6). Hence the performance of the OWSC is essen- 
tially dominated by the exciting torque, as already outlined by [16] . In other 
words, the dynamics of the OWSC is governed by diffractive processes. It is 
then clear that PA and QPA dynamics, governed by radiative processes [5], 
are inadequate to describe the behaviour of large flap-type WECs like the 
OWSC. In figure BJo), the maximum theoretical capture factor Cp ax ( 153]) is 
also plotted for comparison with the actual value Cp (150]) . The two curves are 
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r(s) 



F(s) 



Figure 4: Behaviour of (a) amplitude of rotation in degrees and (b) capture 
factor f )50|) with the period of the incident wave in physical variables, for 
the geometry of table [TJ In (6) the dashed lines represents the maximum 
theoretical capture factor C™ ax (]5B]) . 



very close for short periods. Then C™ ax increases almost linearly in longer 
waves and becomes much larger than Cp. This happens since the OWSC is 
not designed tuned to resonance, i.e. the flap buoyancy torque C and inertia 
/ do not satisfy the resonance condition ( 1521) . As a consequence, v vto > v (see 
149]) and Cp < Cp ax (see[50]). This choice can be justified by considering that, 
if the resonance condition ( )52|) occurred, then the stroke of the flap would 
exceed by far the amplitude of the incident wave, resulting in a situation 
not compatible with the power take-off system [TJ [16] . In conclusion, C^ ax 
(E2J) can be used in short periods for a quick performance assessment of the 
device. Then the analysis must be refined by calculating Cp with the full 
expression f )50|) . In the following section, we shall compare the behaviour of 
the OWSC in the open ocean and in a channel. This will allow us to assess 
the influence of the channel lateral walls on the performance of the system 
in experimental studies conducted in wave tanks. 



4 OWSC in the open ocean vs. OWSC in a 
channel 

Let us compare the behaviour of the OWSC in the open ocean (see §2]) versus 
the behaviour of the same device in a channel of width b' [H see]. Such a 
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1.2 



1.1 



1.0 



0.9 



0.8 




1.2 



1.1 - 



1.0 



0.9 - 



0.8 




Figure 5: Behaviour of (a) exciting torque ratio and (6) capture factor ratio 
between the open ocean model of §|2] and the channel model of Renzi and 
Dias [8] , versus the period of the incident wave. The geometry of the system 
is detailed in tabled] The blockage ratio in the channel model is w = 0.2. 



m 1 


2 


3 4 


5 


T m (s) 9.6 


5.7 


4.5 3.8 


3.4 



Table 2: Resonant periods of the first five sloshing modes for the geometry 
of tabled] in a channel of width h' = 91.6m. 

layout is commonly used in wave tank experiments aiming at reproducing 
the dynamics of the system in the open ocean [T] [TBI see for example]. The 
geometry chosen for comparison is again that of table [U In the channel 
model, the channel width is b' = 91.6m and the relevant blockage ratio 
is w = w'/b' ~ 0.2. The latter corresponds to that used in wave tank 
experiments at Queen's University Belfast on a 1 : 20 scale Oyster 1 WEC 
[TBI d3 f° r details see]. In figure |^a) the ratio \F'JF'\ between the exciting 
torque in the channel \F' C \ and in the ocean \F'\ is plotted against the period 
of the incident wave in physical variables. The effect of the lateral walls in 
the channel configuration is revealed by the spikes occurring at the resonant 
periods of the channel sloshing modes (see table |2]). Such periods correspond 
to the cut-off wavelengths X' m = b'/m, m = 1,2,..., at which the sloshing 
modes turn from propagating to trapped near the device [5] for details see]. 
For the chosen geometry, the influence of the lateral walls results in a 10% 
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4 5 6 7 8 9 10 11 12 13 14 

T'(s) 

Figure 6: Behaviour of the exciting torque ratio between the open ocean 
model of §2] and the channel model of Renzi and Dias j8] for several blockage 
ratios, versus the period of the incident wave. The geometry of the system 
is detailed in table [TJ 

increase of the exciting torque in the channel, at the resonant peaks, with 
respect to the open ocean values (see again figure Eb)- Analogous results 
have been also obtained by Chen [TB] for the horizontal loads on a truncated 
vertical cylinder. Note that a similar behaviour occurs also for the ratio 
Cfc/Cf of the capture factor in the channel Cp c and in the open ocean Cp, 
as depicted in figure Mb)- The increase of the capture factor is of the same 
order of magnitude - about 10% - than the increase of the exciting torque. It 
is then clear that care should be taken when employing the results obtained 
in a wave tank to predict the behaviour of the OWSC in the open ocean, as 
already noted by Renzi and Dias [5] . To quantify the degree of approximation 
of such practice, let us further investigate the effect of changing the width 
of the channel in the wave tank model. Figure [6] shows the behaviour of 
the exciting torque ratio \F' C /F'\ against the period of the incident wave in 
physical variables, for several blockage ratios. When w is small, the channel 
is much larger than the flap and a number of resonant peaks occur (see 
curve w = 0.1 in figure [6]). However, the influence of the lateral walls is 
not appreciable and \F'JF'\ ~ 1. When w increases, the number of resonant 
peaks in the selected period range decreases, while the ratio \F'J F'\ increases 



18 



at peaks. The deviation of the latter ratio from unity is still negligible (about 
10%) for w = 0.2, but is definitely substantial for w > 0.3. Therefore we 
recommend a blockage ratio not larger than 0.2 in wave tank tests aiming at 
reproducing the behaviour of the device in the open ocean. In the following 
section we shall perform a parametric analysis of the OWSC in the open 
ocean with respect to the flap width w'. 

5 Parametric analysis of the OWSC in the 
open ocean 

Here we undertake a parametric analysis of the system to investigate the 
influence of the flap width w' on the behaviour of the OWSC. In figure 
C3a — d) the added inertia torque //, the radiation damping u', the magnitude 
of the exciting torque \F'\ and the capture factor Cf are plotted against the 
period of the incident waves in physical variables, for three different widths 
of the device, namely w' = 12, 18, 26m. The foundation height, the ocean 
depth and the incident wave amplitude are still those of table [TJ Generally, 
increasing the width w' determines all the hydrodynamic characteristics, //, 
v' and \F'\, to increase. The larger the width of the flap, the larger are the 
peaks of said parameters, as shown in figure [7( a — c). A larger flap width, 
in fact, corresponds to a larger surface area, which in turn determines higher 
inertia torque, radiative capacity and exciting torque [UEE see]. Concerning 
the capture factor CV, increasing the width of the flap determines both a 
migration of the peak values towards larger periods and a broadening of the 
curve bandwidth, as shown in figure [T^cf) . The peak values of the capture 
factor curve, however, remain almost constant. This behaviour is different 
from that noticed by Renzi and Dias [8J for an OWSC in a channel. There, 
the wider flaps have also a larger resonant peak, provided the blockage ratio 
is still small [U otherwise the dynamics is quasi-2D and Cf —> 1/2, see]. This 
is due to the selective behaviour of the channel sloshing modes, for which the 
most powerful resonance occurs with large flaps. In the open ocean, on the 
other hand, the radial spreading of energy to infinity inhibits such resonant 
mechanisms and prevents the occurrence of resonant peaks in the capture 
factor curve. Nevertheless, increasing the OWSC width in the open ocean 
has still the beneficial effects of broadening the bandwidth and increasing the 
peak period of curve (see again figure[7]i). In practical applications in random 



19 





20 




16 




12 












8 


=5. 






4 








x W_ 

(a) 



x 10 7 



• w'=26 m 

• w'=18 m 

• w'=12 m 



/ 



/ 



1 



_i I I I l 



10 



4 6 
x 10 6 



10 



12 14 




12 14 



5, 

f 4 

2 




■ (c) 


— • - w'=26 m 

w'=18 m 

-- w'=12 m 








\ 














10 



12 



14 




T (s) 



T (s) 



Figure 7: Behaviour of (a) added inertia torque, (6) radiation damping, (c) 
magintude of the exciting torque and (d) capture factor versus the period of 
the incident wave. The geometry of the system is detailed in table [TJ Three 
different flap widths have been considered. 
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seas, this fact could be exploited by designing the width of the flap so that 
the resulting capture factor bandwidth couples well to the wave spectrum of 
the incident sea. 

6 Conclusions 

A semi-analytical potential flow model is derived for a single OWSC in the 
open ocean. Application of the Green's integral theorem in the fluid domain 
yields a system of hypersingular integral equations for the jump in potential 
between the two sides of the flap. The system is solved via a series expansion 
in terms of the Chebyshev polynomials of the second kind and even order. 
The hydrodynamic parameters of the system are then fully characterised 
and the equation of motion of the flap is solved. Asymptotic expressions of 
the velocity potential are obtained for both the diffracted and radiated wave 
fields. Several relationships are found between the hydrodynamic parameters 
and the amplitude of the radiated and diffracted waves in the far field. Some 
of them are well known in the literature, while others are new and point 
out the distinguishing features of the OWSC with respect to other WECs. 
Analysis for a specific OWSC configuration, inspired by the design of the 
Oyster 1™ WEC, shows that the capture factor of the system is governed 
by the exciting torque acting on the flap. High levels of capture factor can 
be attained, even though the OWSC is not usually tuned to resonance with 
the incident wave field. The behaviour of the OWSC in the open ocean is 
then compared with that in a channel. This allows to quantify the degree of 
accuracy of wave tank experiments aiming at reproducing the behaviour of 
the device in the open ocean. We found that blockage ratios less than 0.2 
guarantee a good agreement between the open-ocean and the channel mod- 
els. Finally, parametric analysis of the OWSC in the open ocean reveals that 
the hydrodynamic characteristics of the system, i.e. exciting torque, added 
inertia torque and radiation damping, all increase by increasing the width 
of the flap w'. The peak capture factor, instead, shows no appreciable vari- 
ations with w' . However, increasing the flap width has the beneficial effect 
of enlarging the bandwidth of the capture factor curve. This phenomenon 
could be exploited in random seas to tune the monochromatic capture factor 
curve to the incident wave spectrum. 

This study was funded by Science Foundation Ireland (SFI) under the re- 
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search project "High-end computational modelling for wave energy systems". 



A Semi-analytical solution 

In this section we find the solution to the radiation and scattering problems 
of §21 by using the Green integral theorem in the fluid domain. The method 
we employ follows the same procedure devised by Renzi and Dias jB] for a flap 
in a channel. Consider the plane potentials (fn' D \ solutions of the system 
02ip - (j2~2"]) and outgoing as r — > oo. Application of the Green integral theorem 
in a large circular region enclosing the flap in the (x, y) plane [HI see] yields 

^R,d) = f 1/2 A^Jfo) G n ,s(x,y;Z, V )\t =0 dr,. (54) 

J -1/2 



In the latter, 



G n (x,y;Z,r,) = (« n y/(x - £) 2 + (y - v) 2 ) (55) 

is the Green function in the (x,y) plane [TD1 see], singular at (x,y) = (£,?]) 
and outgoing as r — > oo, while 

^ D \y) = ^' D) (-o,y) - ^ n R ' D \+o,y) (56) 

is the jump in the plane radiation and diffraction potentials across the plate. 
Following the procedure of [8], apply the boundary condition (l2"2"j) to the 
spatial potentials ( 154"]) . Then expand the Hankel function inside G n (I55I) 
according to [TTJ §8.444], perform the change of variable u = 2rj and define 

(P B ,g B )(«) = (A<p« AvOfa). (57) 



With these manipulations, expression (1221) finally becomes 

h r i ^\ \(v - U )- 2 du 



Qn(u) 

In the latter, which is formally equivalent to (B15) of [5J, the symbol H J 
indicates the Hadamard finite-part integral [8J, [101 see] and 

K>n[f {u)\ = — j— / fiu) f ; du, (59) 

4 J -1 \ v — u \ 
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where v = v (y) = 2y e (—1, 1) and 



RJa) 



Ma) 



2i ( a 
l + -(ln- + 7 



I 

71 



+ g (_iy+i( a / 2 )2i-i 

J=2 



JKJ - 1)! 




(60) 



Ji(x) is the Bessel function of first kind and first order and 7 = 0.577215 . . . 
the Euler constant. Note that K. n fl59|) has a convergent kernel. Hence 
expression (158]) has a singularity only in the finite-part integral of the left- 
hand side. In order to solve ( l58l) in terms of the jump in potentials P n and 
Q n , we seek for solutions of the type 



UlrU=U} (i -" 2)1/2 g{t} 



U 2p (u), 



(61) 



where the a pn and f3 pn are unknown complex constants and the U2 P are the 
Chebyshev polynomials of the second kind and even order 2p, p G N. Sub- 
stituting the series expansion ( 161]) into the singular equation ( 158]) and per- 
forming the algebra yields finally 



E 

p=0 



pn 
pn 



C pn (v ) 



fit 

d„ 



v e (-1,1), 



(62) 



with 



C pn {v ) 



-{p+l)U p {v ) 

Rn(f\v -U\) 



4 ' ^ 



\V - U\ 



du. 



U p {u) 



(63) 



The linear systems of ( 162]) can be now solved numerically by truncating the 
series to a finite number P < 00 and by taking a finite number of evaluation 
points t>o = voj, with j — 0,1, ... P. A fast numerical convergence of (1621) is 
guaranteed when the voj are the zeros of the first-kind Chebyshev polynomials 
\M seel 



(gj + 1)tt . 



0,1,. ..,p, 



(64) 
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for which (K)2"|) reduces to two (P + 1) x (P + 1) truncated algebraic systems 
for each n. The latter can be easily solved to determine the a pn and the (3 pn . 
Finally, some solvability conditions must be satisfied to ensure the uniqueness 
of the solution to f )62|) . For n > 0, d n = (see El]), so that it must be (3 pn = 
to ensure uniqueness of the a pn . Physically, this means that the diffraction 
problem does not admit evanescent modes. For n — 0, instead, the solution 
to (J5"2"|) is unique if 

P P o = <y p od /f . (65) 

The plane radiation and scattering potentials can now be obtained by sub- 
stituting the found jumps in potential (equations [571 and IBTj) into the Green's 
theorem (|54p and by summing up all the vertical eigenmodes according to 

jlZJ- 
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